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The CN0C2 field galaxy redshift survey (hereafter CN0C2) is designed to 
provide measurements of the evolution of galaxies and their clustering over 
the redshift range to 0.7. The sample is spread over four sky patches with 
a total area of about 1.5 square degrees. Here we report preliminary results 
based on two of the sky patches, and the redshift range of 0.15 to 0.55. We find 
that galaxy evolution can be statistically described as nearly pure luminosity 
evolution of early and intermediate SED types, and nearly pure density evolution 
of the late SED types. The correlation of blue galaxies relative to red galaxies 
is similar on large scales but drops by a factor of three on scales less than 
about 0.3h~^ Mpc, approximately the mean scale of virialization. There is a 
clear, but small, 60%, change in clustering with 1.4 mag of luminosity. To 
minimize these population effects in our measurement of clustering evolution, 
we choose galaxies with M^''^ < —20 mag as a population whose members are 
most likely to be conserved with redshift. Remarkably, the evolution of the 
clustered density in proper co-ordinates at r ^ lOh^^ Mpc, pgg oc (1 + z)^, 
is best described as a "de-clustering", oc (1 + 2;)0-6±o.4^ Qj. equivalently, there 
is a weak growth of clustering in co-moving co-ordinates, xq oc (1 + z)~^-^^^-'^. 
This conclusion is supported by the pairwise peculiar velocities which rise 
slightly, but not significantly, into the past. The Cosmic Virial Theorem applied 
to the CN0C2 data gives Qilu/b = 0.11 ± 0.04, where Q is the three point 
correlation parameter and b the bias. Similarly, galaxy groups have a virial 
mass-to-light ratio (evolution corrected) of {M^i^ai / L^^) ~ 215/1L0/M0, or 
VLm = 0.15 ±0.05. 



1. The motivation and design of the survey 

The best test of theories of structure evolution is to observe the evolution of 
both galaxies and their clustering. The Canadian Network for Observational 
Cosmology (CNOC) field galaxy redshift survey is designed to investigate 
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Figure 1. The distribution on the sky, in seconds of arc, of the galaxies with redshifts and 
ruR < 21.7 mag in 19 of the 20 fields in the 0223+00 patch. 



clustering dynamics and its relation to galaxy evolution on scales smaller than 
approximately 20h~^ Mpc over the < z < 0.7 range. To meet these goals 
requires a "CfA-class" survey ( Davis fc Peebles 1983| ) which contains roughly 10^ 
galaxies in 10^ h~^ Mpc^ to allow subsampling and to provide minimal coverage 
of a representative set of clustering environments. For observational convenience 
the survey is distributed over four patches, each nomi nally containing 20 
Multi-Object Spectrograph, MOS, ( [LeFevre et al. 1994 ) fields of approximately 
9' X 8' on the sky. The layout (minus one field) is shown in Figure ||. 

The fields are observed in the UBgRI filt ers, with t he R (Kron-Cousins) 
magnitudes (measured as total magnitudes, |Yee 1991 ) being used to define 
the survey. The photometric sample extends to about mji = 24 mag, with 
comparable depths in the other filters, except for U and I which extend to 
23 mag. The spectroscopic sample is drawn from a "mountain" version of 
the photometric sample with a nominal limit of m/j = 21.5 mag. Each field 
is observed with two spectroscopic masks, the "bright" mask extending to 
mji ~ 20 mag and the "faint" mask extending to the limiting magnitude. Slits 
for additional fainter objects are placed in any otherwise unoccupied areas. 
Together the two masks largely eliminate the problem of slit crowding and 
have the further benefit of increasing the efficiency of the observations. The 
spectra are band limited with a filter extending from 4400-6300A, which gives 
a statistically complete sample over the 0.15 < z < 0.55 range, with emission 
line galaxies visible over < z < 0.7. The average completeness of the resulting 
redshift sample relative to the photometric sample is about 50%. The success 
rate for obtaining redshifts from spectra is about 85% after accounting for the 
fraction of objects expected to be outside our passband. The failures are mainly 
due to poor seeing, poor transparency and objects in the corners of the MOS. 
The redshift distribution for one patch is shown in Figure |2|. The results below 
are derived assuming Ho = 100/ikms^^Mpc"^ and qo = 0.1. 
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Figure 2. The line-of-sight distance versus the transverse distance in the Dec direction from 
the field centre, in proper co-ordinates, in the 0223+00 patch. 



2. Evolution of the Luminosity Function 

The availability of UBgRI photometry for our sample allows us to classify 
CN0C2 galaxies by colour and to compute the luminosity functions (LF) for 
different galaxy populations in a number of different bandpasses. In particular, 
we give details of our LF methods and descriptions of CN0C2 LF evolution 
results in Lin et al. (1998a), and confront our LF, number count, and colour 
distribution data against a variety of galaxy evolution models in Lin et al. 
(1998b). 



In Lin et al. 1998a we calculate LF parameters in the Babi and U bands 



for "early," "intermediate," and "late" CN0C2 galaxies, classified using fits of 



UBRI colours to the galaxy spectral energy distributions (SED's) of Coleman, 
Wu Sz Weedman 198(\ We present a description of the LF evolution using the 



following convenient model, 

M*{z) =M*{0)-Qz 
a{z) = a(0) 
/,(z)=p(0)10°-^^% 

where M* and a are the usual Schechter LF parameters, p is the galaxy number 
density, and P and Q parameterize the rates of number density evolution and 
luminosity evolution, respectively. We plot our Bab LF results in Figure |^ and 
show 2a P-Q error contours in Figure ^. We find that the faint-end slope of 
the LF is steeper for late-type galaxies relative to early-type objects, consistent 



with previous LF studies at both intermediate and low redshifts, (e. g., Lilly 
et g/. 1995| , [Ellis et al. 199^ , |Lin et al. 1996a| ). Moreover, the LF's of the early 



and intermediate populations evolve differently from that of late-type galaxies. 
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Figure 3. Evolution of the Bab luminosity functions (solid curves and points) for early-, inter- 
mediate-, and late- type CNOC2 galaxies in three redshift bins [z increases from top to bottom). 
Also shown are fiducial LF's (dotted curves) from the lowest-redshift bin for each galaxy type. 
Results shown are for go = 0.1. 



Specifically, we find that the LF's of early and intermediate types show primarily 
positive luminosity evolution {Q 1.5) and only modest density evolution 
(P ~ —1), while the late-type LF is best fit by strong positive number density 
evolution (P « 3) and little luminosity evolution {Q !^ 0.5). We also confirm the 
trend seen in previous smaller intermediate-redshift samples that the luminosity 
density of late-type galaxies increases strongly with redshift, but that the 
luminosity density of early- type objects remains relatively constant with z. 
These general conclusions hold for either qq = 0.1 (as in the LF figures shown) 
or qo = 0.5. Specific comparisons against the Canada-France (Lilly et al. 1995) 
and Autofib ( [Ems et al. I996| ) redshift surveys show general agreement among 
our LF evolution results, although there remain some detailed discrepancies 
with respect to the B-selected Autofib survey, which may be due to differences 
in galaxy classification or sample selection methods. 

In p^in et al. 1998a| we also compute SED type distributions, UBgRI number 
counts, and various colour distributions for CN0C2 galaxies. The number 
counts and colour distributions for all R < 21.5 CN0C2 galaxies (not just 
those with redshifts in our completeness range 0.1 < ^; < 0.55 used to compute 
the LF) are well matched once we extrapolate our LF evolution models to 
z ss 0.75, thus providing additional checks on the validity of our LF evolution 
results. In addition, we have verified that various systematic effects, specifically 
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Figure 4. The two sigma confidence contours in P (number density evolution parameter) versus 
Q {M* evolution parameter) for the Bab luminosity functions of early, intermediate, late, and 
early-|-intermediate CN0C2 samples. The intersection of the horizontal and vertical dotted lines 
indicates no-evolution, P — Q — 0. Results shown are for go = 0.1. 

patch-to-patch variations, photometric errors, surface brightness selection, 
redshift incompleteness, and apparent magnitude incompleteness, do not 
significantly affect our results. 

Subsequent papers on galaxy population evolution in CN0C2 will also make 
use of the morphological and spectral information that will become available for 
CN0C2 galaxies once the appropriate data are fully reduced. We are also in 
the process of deriving properly calibrated photometric redshifts, which should 
provide another factor of two increase in useful sample size for R < 21.5 galaxies. 
Future papers will further explore the issue of LF evolution using these even 
larger CN0C2 galaxy data sets. 



3. Two Point Correlations 

On scales where clustering is strongly nonlinear the two point correlation 
function is a dynamically useful statistic, which is now fairly well calibrated 
to initial conditions via n-body simulations (e. g. Jenkins et al. 199^ ). Our 



survey is sufficiently densely sampled, with a velocity accuracy of 100 kms^^ 
or better, that we normally derive all our results from the two dimensional 
correlation function, i{rp,rz), shown in Figure ^, which ratios excess pairs to 
the smooth background at projected separation Vp and redshift separation r^ . 
We estimate i{rp,rz) with the classical estimator DD/DR — 1 (Peebles 1980| ), 
which has the advantage of simplicity and speed. The DD ■ RR/{DR ■ DR) — 1 
and {DD — 2DR + RR) / RR estimators lead to no significant changes of the 
results presented here. 

The real space correlation function can be derived from the projected 

correlation function, Wp{rp) = J ^{Jrp + r'^) dr^- In the limit that the 
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Figure 5. The two dimensional correlation over the 0.15 < z < 0.55 range. The plot is sym- 
metrized about the centre and smoothed with a filter that increases with distance from the 
origin. The apparent flattening of the contours at large distances is not yet statistically signifi- 
cant. 

integral extends to infinity Wp{rp)/rp = r(l/2)r((7 — l)/2)/r(7/2)(ro/rp)'^. The 
immediate complication with this measure of ^(r) is that when applied to the 
data the sum needs to be truncated at some finite r^. The minimum is that 
required to sum over the random pairwise velocities along the line of sight, say 
H{z)rz > 3(Ti2, or approximately Mpc. At larger the sum converges 

roughly as 1 — (ro/r^), that is, relatively slowly. However, increasing the cutoff 
Vz beyond about 50h~^ Mpc increases the noise from large scale structure. The 
results here use cutoff of 20 to 50 Mpc, and fits to the resulting Wp{rp) 
use data with < 10h~^ Mpc. We make no correction for correlation beyond 
the cutoff r^, which is generally less than 10%, even for the unrealistic pure 
power law correlation model. 

(a) Luminosity and Colour Dependence of Correlations 

The auto-correlations of red and blue galaxies (of all luminosities) are shown 
in Figure ^. The sample is split approximately in half at (B — R)q = 1.25 
which divides the galaxies into those with low and high star formation rate. 
The resulting subsamples have nearly identical mean redshifts, 0.35, but the 
red subsample has a mean luminosity of Mji = —20.5 mag as compared to 
the blue M/j = —19.8 mag, so there will be some excess of red over blue as 
a result of luminosity dependent correlation, see Figure |^. The blue galaxy 
auto-correlation has a characteristic scale, roughly 0.3^ Mpc, shortward of 
which they fall well below red galaxies. A similar trend in the Elliptical/Spiral 
ratio is reported in the APM ( [Loveday et al. 1995| ) and CfA-hSSRS2 survey 
(Marzke, private communication). A comparison to the pairwise velocity 
dispersion, approximately 350kms~^, and the properties of groups found in 
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Figure 6. Auto-correlations as a function of colour. The pincushions are for red galaxies, the 
circles are for the blue galaxies. The blue galaxies appear to have a characteristic length for 
correlation change of about 0.3/i~^ Mpc. The error flags are computed from the difference of the 
two fields and greatly exceed the \/iV estimates. 

these data (see below), suggests that the "blue break" at 0.3/i~^ Mpc is the 
radius where field galaxies statistically join a virialized region, inside of which 
star formation is suppressed. The dynamics of virialization in a hierarchical 
cosmology is to tidally remove the outer dark matter and any imbedded baryons. 
The net effect is to eventually starve the galaxy of infalling gas. 

The issue of luminosity dependence of correlations is an important test 
theories of primordial bias associated with density "peaks". The observational 
situation is unclear at the moment. The APM survey finds about a factor of 
t hree increase in correla tion for about a 2 mag change in luminosity. The LCRS 
( iSchectman et al. 19961) has no significant luminosity dependence of clustering 
after allowance for redshift differences, over a directly comparable range of 
luminosities. A mild enhancement, 35% in clustering amplitude, of the clustering 
of high over low luminosity galaxies is seen i n the Perseus-Pisces catalogue 
(IGiovanelli fc Haynes 199l| , |Guzzo et al. 1997| ). 

The dependence of clustering on luminosity is shown in Figure 0. The sample 
is divided at Mr = —20, k corrected and evolution corrected at an approximate 
rate of Mji{0) = Mji — Qz, with Q = 1. The mean luminosities are —19.35 and 
—20.75 mag. The M^, of a Schechter fit to the luminosity function is -20.3 mag. 
The mean redshifts of the two samples are 0.37 and 0.34. The fitted correlation 
lengths are 3.7/i~^ Mpc and 2.65h~^ Mpc both with 7 = 1.7, or a ratio of 
clustering amplitudes of high to low luminosity of 1.67. These correlations are 
significantly larger than estimated from the small sky area, somewhat lower 
lumin osity, samples available previously (p^eFevre et al. 1996 , Shepherd et al. 
199^ , ICarlberg, Cowie, Songaila fc Hu 199'^ . 



Within the CN0C2 data the change of correlation amplitude with luminosity 
is comparable in size with what is expected on the basis of galaxies of increasing 
luminosities being correlated with dark matter halos of increasing mass. That 
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Figure 7. Auto-correlations as a function of luminosity. The pincushions are for A/^'' 
galaxies and the circles for those between —18.5 and —20 mag. 



< -20 



is, very approximately, ^ cx A(M) ^ oc M^"'~^^^^^, where A(M) is the mass field 
variance for spheres containing mass M, which for a perturbation spectrum, 
P{k) cx /c", is A cx M"*^""*"^)/^. For n ~ — 2 expected from CDM-like spectra 
on galaxy scales, we then expect a correlation amplitude ratio of 1.58 for the 
factor of four difference in luminosity. That is, we attribute the luminosity 
dependence of the correlation as reflecting an underlying primordial difference 
in the correlations, roughly as expected for galaxy formation in dark matter 
potential wells. 



4. Evolution of the Two Point Correlation Function 

At low redshift the evolution of galaxy correlations can be adequately 
described with S,{r,z) = (ro/r)'^(l + z)~^'^~^''\ where the lengths are measured 
in proper units. This double power law model does not allow any variation 
of the correlation slope, 7, with redshift. The model might seem theoretically 
naive, but it is usually a better description of the data than any available 



non-linear realization of a range of CDM models (Colin, Carlberg Sz Couchman 



1997 , Jenkins et al. 1997). A t low redshift 7 = 1.8 is the canonical empirical 



value ([Davis fc Peebles 1983| ). As redshift incre ases there are three very gener al 



possibilities as, verified in n-body simulations ( Colin, Carlberg Sz Couchman 
|T99^ ): 

• e ~ for particle clustering if r2j\/ low, 

• e ~ 1 for particle clustering if — 1, and, 

• e ~ — 1 for dark matter halo clustering, with weak dependence. 
The entire set of n-body particles defines the evolution of the mass field. Halos 
with densities of 200pc evolve in numbers as described by the Press-Schechter 



model ( Press &: Schechter 1974 ). However these densities are too low to work 



as a one-to-one association of galaxies inside groups and clusters, although 
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Figure 8. The BC96 model estimated of the specific star formation rate (multiplied with 15Gyr, 
assuming a duty cycle of unity) versus the model mass for the CN0C2 galaxies. Galaxies of low 
mass and low star formation rate will be below the m_R limit of our survey. 

they may be appropriate for less clustered field galaxies. Halo cores at densities 
well abo ve the virializat i on density, 10^ — W ^Pc, remain distinct even in large 



clusters ( Carlberg 1994 , Ghigna et al. 199^ ) and are the best candidates for 
identification with galaxies in the context of purely collisionless large scale 
structure simulations. 

A primary issue in studying the evolution of correlations is to be able to 
identify the same population at different redshifts, since both luminosity and 
colour dependence of clustering can mimic or mask the desired effect. Although 
not yet fully implemented, fitting the photometric SED to a model gives an 
indicative mass-to-light ratio which can then be used estimate galaxy masses, as 
is shown in Figure |. We use the Bruzual and Chariot (1996) SED r = 0.5, 1, 
1.5, 2, 4 and 20 Gyr models and ages from 1 to 19 Gyr in steps of 2 Gyr. The 
very blue galaxies require 1 Gyr old models, which is the cause of the gap in 
the indicative star formation rate of Figure ^. The fitted models give a stellar 
M/L and a specific star formation rate, which is multiplied by 10 Gyr to give 
an indication of the significance of star formation over a Hubble time, assuming 
that the duty cycle for star formation is 100%. There are two indicative results 
from Figure ^ which will motivate our sample choice for correlation evolution 
estimation. First, the galaxy mass function shows no significant evolution over 
our redshift range, although a 30% or so variation would be within the limits. 
Second, there is a clear confirmation of the strong statistical correlation between 
star formation and the stellar mass of galaxies. 

High luminosity galaxies, M^'^ < —20, are the closer than lower luminosity 
galaxies to being a certifiable mass invariant population, although the bluest 
members of this set likely have significantly lower stellar masses which may 
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Figure 9. Estimates of ro and 7 for high luminosity galaxies from the LCRS and CNOC2 samples. 
The mean redshift, from top to bottom is 0.080, 0.135, 0.28, and 0.43. Contours are 68, 90 and 
99% confidence levels. 

not be well conserved over our redshift range, if their star formation is steady 
in time. This population has the considerable advantage that an identically 
defined sample can be found in the LCRS. The high luminosity galaxies with 
< —20 comprise a volume limited sample within the CN0C2 data, and 
approximate one within the LCRS sample. 

Two point correlation parameters, rg and e, for a power-law correlation 
function of the high luminosity galaxies are shown for various redshifts 



Figure 0. We find that tq = 5.15 ± 0.15/1"^ Mpc and e = -0.6 ± 0.4. 

These data strongly exclude clustering evolution that declines as rapidly as 
e = 1. One could erroneously infer such a rapid decline if one used a sample in 
which galaxies at higher redshifts are intrinsically less luminous or on the average 
bluer than those nearby, as is the case for the general population as a function 
of redshift. The rate of correlation evolution falls between the values expected 
for low Qj\,f particles (non- merging) and dark matter halos (which do merge). 
Our result can be re-stated in terms of the evolution of the correlation length in 
comoving co-ordinates, xq = ro(l -|- z)^(^+'^^'>')/'>'), as xq cx {1 + z)~^-^^^-'^. 

It should be borne in mind that our estimate of correlation evolution is a 
preliminary result and that the errors will be reduced with the full sample. 
An important consideration is that the luminosity cut that defines our sample 
mixes together galaxies of a fairly wide range of masses and anti-correlated star 
formation rates. It seems likely that low and high star formation rate galaxies 
have different correlation histories, which are mixed in the current approach. 



The redshift space distortions in (,{rp,rz) reflect the dynamics of clustering. 
The random velocities elongate contours of £,{rp,rz) at small and infall 
velocities squash the contours at large r^. These distortions depend on 0,, and 
the biasing of galaxies with respect to the mass field. The data here have a 
local velocity accuracy of better than 100 kms~^, as explicitly demonstrated 
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Figure 10. Estimates of ro at z = and e for high luminosity galaxies from the LCRS and 
CNOC2 samples. The result falls between the e = (fixed physical clustering) and e = —1.2 of 
clustering fixed in co-moving co-ordinates. 
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Figure 11. The versus model pairwise peculiar velocity, ai2 within the LCRS sample (solid 
line) at a mean redshift 0.10 and the CN0C2 sample (dashed line) at a mean redshift of 0.36. 



using redundant spectra taken through different spectrograph slit masks. This 
is comparable to the velocity accuracy of surveys at low redshift. The CN0C2 
survey is designed to concentrate on scales less than lOh^^ Mpc so will not 
provide strong limits on the infall velocities. 

The pairwise peculiar velocities are derived from a model for £,{rp, r^) following 
the procedures of Croft, Dalton & Efstathiou (1998). Once the rp and 7 are 
derived from the velocity independent correlations, then we set VLm = 0.2 and 
compute as a function of ai2- In Figure 11 we show the reduced versus 
the pairwise peculiar velocities for the LCRS sample at a mean redshift of 0.10 
and the CN0C2 sample at mean redshift 0.36. The minimum of rises slightly 
with redshift, although it is consistent with no change with redshift. The e 
model for clustering predicts with the Cosmic Virial Theorem (see below) that 
(Ti2(z) oc (1 -|- z)"*^/^, provided that the bias is not changing with redshift. We 
conclude that the peculiar velocities evolve in accord with that predicted from 
the correlation function alone. This is strong evidence that biasing of galaxies 
with respect to dark matter is not a large effect at these redshifts. 

The peculiar velocities also show a strong population dependence, as shown 



in Figure 
whereas t 



j. The red galaxies have a pairwise dispersion of about 350 km s 

ae blue galaxies have a measured dispersion of about 200kms~^. 
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Figure 12. The versus model pairwise peculiar velocity, rju for the red (solid line) and blue 
galaxies (dashed line) within the CNOC2 sample. There is between 70 and 100 kms~^ added 
in quadrature with the true velocities. 

Removing the velocity errors in quadrature would reduce the pairwise velocities 
of blue galaxies to about 150kms~^, which is a remarkably cold population. One 
can speculate that in as much as accretion is necessary to promote ongoing star 
formation, then star formation should be suppressed in regions of strong tidal 
fields (the groups) and promoted in moderately dense regions of low velocity 
dispersion. A plot of the estimated star formation rate versus local phase density 
indicates a weak, but suggestive correlation. 

The cosmic virial theorem (CVT) estimates the value of in the field 
( Pavis &: Peebles 1983 ) and is an important complement to studies of clusters 



as Q,M indicators. For a power law ^(r), the CVT reads 

Mm <jh 4(7-l)(2-7)(4-7) 



Q- 



b (1 + z)^r^r^-^H^ 3 J(7) 



( pavis fc Peebles 1983| ) where Q is the three point correlation parameter, b the 
linear bi as of galaxy c lustering relative to mass clustering, and J{'j) is 4.14 for 



7 = 1.7 ([Peebles 19801) . For au = 350 ± 50, ro = 3.2 ± 0.2 (the mean val ues ove r 
the 0.15 < z < 55 range) we fin d QQM/b = 0.11 ± 0.04. For Q = 0.7 (|Di^ 



Peebles 198^ , [Bean et al T983|) this indicat es 0.15 ± 0.06. This is 



in good accord with the results from clusters ([Carlberg et al. 1996| , Carlberg, 



Yee & Ellingson 1997) and the galaxy groups within this survey, neither of 
which is significantly affected by bias. Because our sample finds both a uniform 
luminosity dependence of clustering and a scale dependent colour effect, all 
galaxies cannot have 6 = 1. However, it is clear that on lOh^^ Mpc scales 
and smaller, that values of 6 > 2 or 6 < 1/2 are very unlikely. Moreover, this 
is further evidence that the matter density of the universe is ~ 0.2 and 
possibly somewhat lower. 



6. Galaxy Groups 

Galaxy groups are the dense, quasi- virialized, regions within the field. As 
virialized regions they offer several interesting possibilities: 

• they can be used to estimate il.M via Oort's method, 

• the mean mass and luminosity profiles can be independently measured, 

and. 
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Figure 13. The number of groups found as a function of their hne of sight (rest frame) velocity 
dispersion. The groups are selected from a volume of about 2 x 10^ h~"^ Mpc''. Below 100 kms~^ 
there are few groups due to our velocity measurement error. 



• the relation between clustering and galaxy properties can be investigated. 
The highly nontrivial problem is to identify groups using redshift space data. We 
adopt a simple algorithm. We ratio the number of (selection function weighted) 
galaxies within a "window" to those in a random sample, as for the correlation 
function, to measure the redshift space density at the location of every galaxy. 
The window is generally 0.5h~^ Mpc in radius and ±600 kms~^ in the redshift 
space direction. We then select from this list the galaxy with the highest 
overdensity and join to it all galaxies within a second, slightly larger window, 
1200 kms~^ and 0.8 h^^ Mpc. The sole output from this operation is an estimate 
of the location of the centre of the group in RA, Dec and z. The resulting 
centre has quite a weak dependence on the details of the windows provided 
they are not so small that groups are overlooked, or, that they are so big that 
huge fluffy regions are selected as groups. The velocity cutoff is about three 
times the pairwise velocity dispersion. The cutoff in projected radius is 20-30% 
of the correlation length, approximately the radius at which the redshift space 
correlation function, ^(s), bends due to random velocities. Our estimated centre 
is then given to an algorithm which estimates the velocity dispersion, mass and 
total luminosity from galaxies within some aperture, usually OAh~^ Mpc and 
900 kms~^. These values are not particularly crucial, since neither the projected 
velocity dispersion nor the inferred M/L has a measurable radial gradient. 

The groups are shown in Figure 13 as a function of their line of sight velocity 
dispersion. Although these groups should have only a small population bias to 
estimate the global M/L, but they need to be calibrated against n-body data 
if their cosmological number densities are to be compared to Press-Schechter 
predictions. The comforting news from this figure is that there is no evidence for 
many high velocity dispersion groups which are likely to be bogus. If anything, 
it is likely that the ~ 10% of galaxies assigned to groups here means that we are 
missing some genuine groups. 

The virial mass-to-light ratio of the groups with six or more members is shown 
as a function of velocity dispersion in Figure 14. The groups as a whole have a 
mean mass-to-light of about 215/iM0/L0, and the groups as individuals appear 
to be nearly consistent with having the same value. On the other hand, there 
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Figure 14. The group mass-to-light ratio as a function of velocity dispersion for groups with six 
or more members. The drop in M/L at smaU velocity dispersion may be a result of correlated 
errors in both axes, however it could have an origin in dynamical friction. 
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Figure 15. Group luminosity plotted against mass. These quantities are not intrinsically corre- 
lated by the measurement. The solid line is the average group MjL = 215/iM0/L0 and the 
dashed line is for the population adjusted rich cluster MjL = 300/i / . The more massive 
groups appear to be consistent with clusters, but the smaller ones appear to have lower M/L 
values. 

is an indication that lower velocity dispersion groups have lower MjL values. 
The two axes in Figure |lj are statistically correlated, therefore we plot the 
mathematically independent quantities of group luminosity versus virial mass in 
Figure |l^. The more massive (or equivalently, higher velocity dispersion) groups 
have MjL values consistent with the mean value of SOO/iMq/Lq for CNOCl 
clusters, after adjusting 0.12 mag to account for population differences and the 
slightly different luminosity scales. An indication that lower luminosity groups 
are "under-massive" is still apparent in Figure This is perhaps not too 
surprising because there must be a scale at which the MjL characteristic of the 
mean field, SOO/iMq/Lq or so, drops into the "linear rise" regime characteristic 
of individual isolated galaxies ( [Zaritsky &: White 19941 ) . 
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Figure 16. The mean surface density of the groups (arbitrary units) versus radius. The ratio 
of red (dashed) to blue (dotted) galaxies increases towards the centre. The solid line is for all 
galaxies. 

7. Internal Properties of Groups 

The cosmic virial tiieorem uses pairwise statistics, so its il.M is dependent on 
a knowledge of the bias of the tracer galaxies relative to the mass field. On the 
other hand, if the single particle velocity distribution can be found it allows the 
mass field to be derived independently of the distribution of the tracer galaxies. 
Galaxy groups provided that opportunity, since the velocities can be measured 
with respect to the group mean velocity. This means that we can put limits on 
the bias of galaxies with respect to the mass field on scales of about 0.5h~^ Mpc 
and less. 

We use the groups whose local centers indicated velocity dispersions more 
than 150 kms~^ and had six or more members. The following analysis uses 
the group mean center (where there is usually no galaxy so the density reaches 
a plateau) and includes all galaxies within 600 kms~^ in the frame of the 
group. Jeans Equation allows the mass field, M{R) = J p{r) dV, to derived from 
a tracer population, i^(r), which must be effectively in equilibrium, but not 
necessarily distributed like the mass. 



M(r) 



din a?. dlnu 

H 

d In r d In r 
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One complication is the velocity anisotropy parameter, /? = 1 — ag/a'^, which 
has a large effect near the center but this uncertainty diminishes to about 20% 
near the virial radius (Carlberg, Yee & Ellingson 1997). 

The average group surface density profile, T,{R), Figure |l^, is based on about 
200 galaxies. The mean projected slope between logi? = —1 and logR = —0.5 
of all the galaxies is almost exactly dlog(S)/(ilog(r) = —1, which deprojects 
to u (X of the isothermal sphere. The red galaxies are more concentrated 
than the blue galaxies in these groups, with small differences in the luminosity 
function from the field. The projected velocity dispersion profile. Figure 17, 
of both the full sample is consistent with no radial variation in the mass to 
total light ratio. However the blue galaxies have significantly higher velocity 
dispersions than the red galaxies with an indication of a drop with radius. 
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Figure 17. The projected velocity dispersion versus radius in the groups. Blue galaxies (dotted) 
have systematically higher velocity dispersions than the red galaxies (dashed) , which are similar 
to the full sample (solid). 

This is quite distinct from the pairwise velocities, in which blue galaxies have 
significantly lower velocities than red galaxies. This helps one to understand 
why the virial masses given by blue galaxies are about 50% higher than the 
full sample. Red galaxies give virial masses that are 98% of the full sample 
result. We take these surface density and velocity dispersion profiles as tentative 
evidence that the full galaxy distribution and the mass distribution in groups 
are similar. Within groups, blue galaxies are "anti-biased" with respect to the 
red galaxies and likely the mass as well. 

Although these results are both preliminary and clearly subject to statistical 
uncertainties, they would support the notion that these groups are in many ways 
dynamically similar to rich clusters although they do not have such extreme 
population differences as rich clusters and the field. These issues are a major 
area of investigation in the CN0C2 survey. 



8. Ultra-Large Scale Power 

Measuring the power spectrum on scales of lOO-500/i^^ Mpc, roughly the scale 
of the "first Doppler peak" in the CMB fluctuation spectrum, requires a survey 
that covers an appreciable fraction of the visible universe. Pencil beam surveys 
do just that, but only in one dimension which leads to the problem of aliasing 
in of shorter scale power into the derived ID power spectrum, Pi{k). In view of 



previous results (Broadhurst et al. 1990) and the fundamental interest in this 



statistic we measure the one dimensional power spectrum from the CN0C2 data 
and compare it to various simple models. 

The one dimensional power spectrum Pi{k) = \N~'^Y^f_,eyiY>[ikxi\\^ 
(unweighted) is shown as the upper irregular line in Figure |l^ At very large kz 
this tends to 1/A^ where N ~ 1500 in each of our patches. The observed ID 
spectrum requires some care in interpretation since the ID power at a large 



wavenumber is an integral over the entire 3D spectrum (Kaiser & Peacock 1991) 
whose analysis we follow here. The ID power spectrum of the smooth n{z) is 
the W{kz) window function and is shown as the lower irregular line. Clearly 

the power at kz < 0.01h~^ Mpc is completely dominated by the smooth 
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Figure 18. The one dimensional power spectrum and related quantities. The lower irregular line 
is the power spectrum of the smooth n{z) distribution. The upper irregular hne is the power 
spectrum of the observed redshifts averaged over the two fields. The upper smooth line is for 
{'''o/rY''^ , with ro = 4.2h~^ Mpc (co-moving) which is modified with a Gaussian peculiar velocity 
distribution of 300 kms~^ and then the power correlation truncated at 50 and 100 Mpc. 

distribution of the data. At shorter wavelengths, A < 600/i^^ Mpc, the observed 
power spectrum is a result of real structure. 

The null hypothesis is that the ID power spectrum is simply the power 
spectrum of the real space correlation function, (ro/r)"^, as projected into 
these observations. The 3D power spectrum is P{k) = 47rA;~^(A;ro)T(2 — 
7) sin [(2 — 7)7r/2]. To derive the window function on the sky we do a Monte 
Carlo evaluation of W{kx,ky) by randomly populating the sky patches and 
transforming. To evaluate the predicted ID spectrum we approximate the z 
window function as 2tt/Lz, and integrate the product of the sky window function 
with P{k). The results are displayed in Figure |l^, where we have added a 
constant shot noise of 1/A'^. The unaltered spectrum does a very poor job of 
describing the observed spectrum. The two necessary modifications are to add a 
random peculiar velocity and a long wavelength decline in power. The random 
velocities are modeled using a Gaussian (clearly not a very good description of 
these data) at 300 kms~^. Any CDM type spectrum rolls over in the 30-100 
Mpc range. We simply truncate the input power spectrum at 50 and 100 
Mpc to produce the plots displayed. We can conclude that the universe is 
not a fractal with a fixed clustering dimension to all scales, and is otherwise 
consistent (at the current statistic al level) with the known large scale clustering 
power spectrum ( [Lin et al. 1996b ). 

9. Preliminary Conclusions 

The CN0C2 sample spans a range of redshifts which has dramatic galaxy 
evolution but in which detailed studies are readily performed. Our preliminary 
conclusions are discussed below. 

Galaxy evolution is confirmed to be strongly differential with the blue, 
low mass galaxies changing significantly in their visible numbers. We find 
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Figure 19. The evolution of the proper density of galaxies around galaxies, with the Ad'^" LCRS 
and CNOC2 sample at low redshift and the Giavalisco et al. (1998) inversion at high redshift. 



luminosity functions similar in their overall features to previous studies ( p^illy 
et al. 1995, pllis et al. 1996 ) but with smaller statistical errors. The evolution 
of the early and intermediate SED types is well described as a pure luminosity 
evolution. In contrast, the late type SED population is described as nearly pure 
density evolution. It should be noted that the density evolution could be at least 
partly in the form of varying duty cycle of star formation. 

The correlation of galaxies above M^, + 0.3 mag is about 60% higher than 
those below, with a mean sample luminosity difference of 1.4 mag. This is 
consistent with the expected intrinsic correlation differences of dark matter halos 
of different mass scales. Having a cleaner mass cut (from the SED model fits) 
and more data to trace the evolution of this effect will allow a more confident 
conclusion. 

Blue galaxies are about 50% less correlated than red galaxies on scales 
beyond about 0.3/i~^ Mpc. On shorter scales blue galaxies have a break in their 
correlation function and become about a factor of three less correlated than red 
galaxies. On large scales the difference is about what one would expect from 
the smaller average masses of the blue galaxies. Because the ^(r) break occurs 
at about the average virialization scale it suggests that the blue-red difference 
is an active environmental effect. Dynamically virialization occurs where tidal 
stripping begins to remove the dark halos and their contents from infalling 
galaxies. 

A primary result is that we find that the correlation evolution of galaxies 
brighter than Mr = —20, where = —20.3 mag, are well described with 
(^(r) = (ro/r)'^(l + z)"^"*^, where tq = 5.15 =b 0.15/i~^ Mpc (proper co-ordinates), 
7 = 1.77 ± 0.05, and e = —0.6 ± 0.4. That is, the physical density of galaxies 
around galaxies is declining slightly with redshift. Figure |l^ shows that the 
extrapolation of our rate of clustering evolution passes smoothly through the 
inversion of the angular correla tion of th e high luminosity z ^ 3 "Steidel 
objects" ( Giavalisco et al. 1998 , but note Adelberger et al. 1998| ). We note that 
the association of our population of objects with theirs is not secure. 

The measured e is completely incompatible with the e ~ 1 expected for the 
mass field in an ^Im = 1 universe. The result is about half way between mass 
clustering in a low density universe and the clustering evolution of the relatively 
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high density cores of dark matter halos ( Pohn, Carlberg &: Couchman 1997 ). 



In general an e < imphes that bias is changing with redshift, although in a 
low Qm model this is a very slow function of redshift. An interpretation of the 
slow drop in proper density, within the context of n-body simulations, is that it 
is due to merging of galaxies in high density regions. Taking the result at face 
value, this means that about 35 ± 20% of all the high luminosity galaxies will 
have merged since redshift one. If the time scale for merging is 1 Gyr, then the 
merging fraction will be about 4 it 3%, on average, close to the estimates of 
Patton et al. (1997). 

The pairwise velocity dispersions of the luminous galaxies do not detectably 
evolve between z = 0.05 and z = 0.4, remaining constant at about 350 kms~^ (for 
our chosen estimator). The Cosmic Virial Theorem gives Q^m/^ = 0.11 ± 0.04, 
where Q is the three point correlation parameter. The three point parameter 
will be estimated from these data, with the expectation that Q — 1 from low 
redshift results. 

Galaxy groups with line-of-sight velocity dispersions between 100 and 400 
kms~^ are found, with numbers declining with increasing cJi, very roughly 
as expected from Press-Schechter theory. The virial M/L values may drop at 
small o"! (perhaps due to dynamical inspiral) with the higher velocity dispersion 
objects having M/L values consistent with the population adjusted values of 
rich clusters. The mean M/L is about 35% below rich clusters. The groups are 
consistent with the total galaxy light being distributed like the mass. The groups 
have a weak population gradient, becoming redder from the "edge" of the mean 
group toward its centre. We tentatively associate the break in the £,{r) of blue 
galaxies with this behaviour of groups. The groups, averaged together, indicate 
= 0.15 ±0.05. 

The one-dimensional power spectrum of our pencil beams requires that the 
the observed power law correlation at scales of less than 10/i~^ Mpc must have a 
turnover at about lOO/i^^ Mpc, which is of course expected in CDM-like theories 
and indicated by observations of large scale clustering at low redshift. 

We set out to test theories of galaxy and structure evolution by measuring the 
evolution. In spite of the dramatic changes in the numbers of lower luminosity, 
blue galaxies, the systems of higher luminosity that dominate the stellar mass of 
the universe appear to be nearly a static population, having only a small change 
in the density of similar galaxies around them and little star formation over our 
redshift range. 

These preliminary analyses are based on 1/2 of the eventual data for which we 
complete observations in May 1998. The results are likely to change somewhat 
(especially as other estimators are implemented and systematic errors detected 
and eliminated). Perhaps the primary conclusion to be drawn is that the CN0C2 
dataset will be a rich sample for the study of cosmology, large scale structure 
and galaxy evolution. 
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